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ABSTRACT 

We propose a new numerical procedure to simulate a single dark halo of any size and 
mass in a hierarchical framework coupling the extended Press- Schechter formalism 
(EPSF) to N-body simulations. The procedure consists of assigning cosmological ini- 
tial conditions to the particles of a single halo with a EPSF technique and following 
only the dynamical evolution using a serial N-body code. The computational box is 
fixed with a side of 0.5/i^^ Mpc. This allows to simulate galaxy cluster halos using 
appropriate scaling relations, to ensure savings in computing time and code speed. 
The code can describe the properties of halos composed of coUisionless or coUisional 
dark matter. For coUisionless Cold Dark Matter (CDM) particles the NFW profile 
is reproduced for galactic halos as well as galaxy cluster halos. Using this numerical 
technique we study some characteristics of halos assumed to be isolated or placed in 
a cosmological context in presence of weak self-interacting dark matter: the soft core 
formation and the core collapse. The self-interacting dark matter cross section per 
unit mass is assumed to be inversely proportional to the particle collision velocity: 
a/rrix oc 1/v. 

Key words: cosmology - dark matter-galaxies: haloes. 
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1 INTRODUCTION 

A detailed understanding of structure formation is one of 
the central goals of contemporary astrophysics and cosmol- 
ogy. In the popular hierarchical clustering framework (White 
& Reese 1978) luminous galaxies form by gas cooling and 
condensing within dark matter halos. These halos merge to 
build larger structures, hierarchically. In the CDM scenario 
systematic studies of halo density profiles for a wide range of 
halo masses were derived by Navarro, Frenk & White (1997; 
hereafter NFW) who argued that the analytical profile of 
the form; p(r) — ps{r /rs)~^ {1 + r/rs)~'^ provides a good 
description of halo profiles for all halo masses, where rs is 
the scale radius which corresponds to the scale at which the 
slope of the profile is —2. 

While predictions of the CDM models have successfully 
accounted for many observations at large scale scales, in- 
creasing interest in testing the predictions at subgalactic 
and galactic scales has grown in the last few years. Inter- 
est began from indications that observed HI rotation curves 
in the central regions of dark matter dominated dwarf and 
Low Surface Brightness (LSB) galaxies are at odds with pre- 
dictions from hierarchical models (Flores & Primack 1994; 



Moore 1994; Burkert 1995). The models predict circular ve- 
locities that increase too rapidly with growing radius com- 
pared to the observed profiles of rotation curves. This fact 
implies a failure of the shape of the predicted halo den- 
sity profile. Recent high-resolution N-body simulations have 
shown that, as the numerical resolution is increased, the in- 
ner profiles is p oc r~^'^ , even steeper than the NFW profile 
(Moore et al. 1999; Fukushige & Makino 2001), exacerbat- 
ing even more the discrepancy with the observations. On 
the other hand, Swaters, Madore & Trewella (2000), provid- 
ing Ha rotation curves for some LSB galaxies, challenged 
the existence of soft cores for these galaxies. They pointed 
out that the HI rotation curves are affected by poor spatial 
resolution, causing a smoothing of the curves in the inner 
regions (usually named beam smearing). For a sample of 
late- type dwarf galaxies, van den Bosch & Swaters (2001) 
show that beam smearing does not allow for distinguish- 
ing between the existence or lack of soft cores in late-type 
dwarf galaxies. However, high-resolution Ha rotation curves 
of LSB galaxies of de Blok, McGaugh & Rubin (2001) and 
Marchesini et al. (2002) are in favour of core-dominated ha- 
los. Similar results are derived by Salucci & Burkert (2000) . 
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While the question of the existence or lack of shallow cores in 
galaxies continues to be debated, some authors have revealed 
the presence of soft cores at the centre of some clusters of 
galaxies from strong lensing observations in C10024+1654 
(Tyson, Kochansky & DoU'Antonio 1998) and from X-ray 
data in Abell 1795 (Ettori & Fabian 2002). Intriguely, clus- 
ters showing a soft core do not have a prominent central 
cD galaxy. In particular Abell 1795 has a central cD galaxy 
moving far from the centre (Ettori & Fabian 2002). On the 
other hand, Chandra X-ray data show the lack of any core 
at the centre of HydraA (David et al. 2001) and a very small 
core in EMSS 1358+6245 (Arabadjis et al. 2002); both clus- 
ters have a central cD galaxy. In relaxed clusters, a dominant 
central galaxy could play a role making steeper the primor- 
dial halo density profile. Accepting the presence of soft cores 
advocated in Abell 1795 and C10024-I-1654, Firmani and co- 
workers (2001) have correlated the halo scales from galax- 
ies to galaxy clusters, integrating the available information. 
Surprisingly, they find that, in this case, the central den- 
sity of dark halos is independent of the halo mass with a 
value of 0.05/i^ Mq pc~^, and that the core radius increases 
proportionally with the halo maximum rotation velocity. 

The presence or the lack of soft cores at the centre of 
galaxy clusters has interesting consequences for the nature 
of the dark matter. The lack of any core in clusters makes 
attractive several scenarios in the debate over the nature of 
the dark matter. In warm dark matter models, the maximum 
phase space density of the particles defined as /max = Po/v'^, 
with po the halo central density and v the halo velocity dis- 
persion, has a finite value and is preserved due to the Liou- 
ville theorem, implying an increase of the halo central den- 
sity when the halo mass increases: po oc (Sellwood 2000; 
Hogan & Dalcanton 2000). In this scenario very massive 
halos of galaxy clusters are predicted to have high central 
densities. 

A different nature for the dark matter was proposed by 
Spergel & Steinhardt (2000) in order to overcome the soft 
core question: if the dark matter is self-interacting, heat 
transfer towards the central regions triggers a thermalizar 
tion process in the dark halos avoiding the formation of a 
cuspy profile. The effects of weak self-interacting dark mat- 
ter were investigated using N-body simulations on isolated 
halos (Burkert 2000; Kochaneck & White 2000) and on CDM 
halos (Yoshida et al. 2000; Dave et al. 2001). These simu- 
lations were performed assuming a cross section indepen- 
dent on the relative particle velocity. Ostriker (2000) and 
Hennawi and Ostriker (2002) have shown a possible incon- 
sistency of the coUisional scenario: indeed the model would 
cause an exorbitant grow on supermassive black holes that 
imposes a very strict upper limit on the collision cross sec- 
tion. This conclusion derives from their assumption that 
when black holes seeds have been formed (at z < 20) the in- 
nermost dark halo shown a NEW density profile. This is not 
necessarily true, at z < 20 collisions have lowered the cen- 
tral density to values at which the growth of the black holes 
seeds becomes negligible. The other upper limits pointed out 
by the same authors and concerning the gravothermal catas- 
trophe and the galactic halo evaporation in clusters are con- 
sistent with the cross section used in this work. On the other 
hand, several authors have pointed out some limitations of 
the coUisional dark matter models. Assuming a dark mat- 
ter cross section independent of the particle relative veloc- 



ity, Miralda-Escudc (2002) shown that this scenario predicts 
cluster cores which are too large and round to be consistent 
with gravitational lensing data. 

This work is focussed on the soft core question and ex- 
plores the possibility that the cross section for the dark mat- 
ter interaction decreases with velocity a/rux oc 1/v as sug- 
gested by Firmani et al. 2001 and Wyithe et al. 2001. The 
focus here is an exploration of halo core properties using N- 
body techniques in a very weak self-interacting regime. For 
haloes within a hierarchical context, the cross section value 
we propose rules out ranges in which the evaporation prob- 
lem and the core collapse could be significant. We show that 
if further observational data can confirm evidence of shallow 
cores in galaxy cluster halos, then self-interacting dark mat- 
ter with a cross section inversely proportional to the halo 
dispersion velocity is capable of predicting the existence of 
soft cores in dwarf galaxies as well as in galaxy clusters. This 
paper is organized as follows: in section 2 the numerical tech- 
nique is described, in section 3 for coUisionless dark matter 
particles, we compare the cosmological halo density profiles 
obtained with our technique to the NFW models for differ- 
ent size halos. Section 4 introduces the Monte Carlo method 
for weakly interacting particle systems implemented in the 
code. In sections 5 and 6, the halo density profiles are in- 
vestigated in a self-interacting dark matter scenario for two 
different cases: isolated halos and cosmological halos. 



2 METHODOLOGY 

In all numerical simulations, approximations and compro- 
mises are necessary if the calculation is to be completed 
within a reasonable amount of time. An important choice to 
make is the effective resolution of the simulation: increasing 
the spatial resolution requires increasing both the number 
of particles to prevent two-body effects and the number of 
timesteps to follow the evolution of smaller structures. For a 
given simulated volume of universe, the number of particles 
determines the mass resolution. In cosmological simulations 
involving only coUisionless dark matter the choice of particle 
number is often driven by computer memory limitations. At 
present large cosmological N-body simulations have reached 
the stage where detailed structural properties of many dark 
matter halos can be resolved simultaneously making use of 
codes running in parallel and of a very large and expensive 
number of particles. We use for our simulations the adap- 
tive, particle-particle, particle-mesh P^M-SPH code, which 
is basically the serial publicly released version of HYDRA 
described in detail by Couchman, Thomas & Pearce (1995). 
The publicly available serial version of the code does not 
allow for cosmological simulations of large volumes of uni- 
verse and the resolution of smaller structures. With the aim 
of overcoming this difficulty without using supercomputers, 
we have introduced relevant modifications to the code in or- 
der to satisfy the following requirements: i) the code has to 
be reasonably fast, serial, and running on a single processor 
workstation; ii) for each run the code must follow the dy- 
namical evolution of a single isolated or cosmological halo in 
a limited volume; lii) the code has to be capable of describ- 
ing the properties of coUisional or coUisionless dark matter 
halos of different sizes and masses from galactic to galaxy 
cluster scales with an arbitrary number of particles. 
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The main change to the code concerns the cosmolog- 
ical initial conditions. Generally cosmological simulations 
with N-body grid codes begin by fixing a power spectrum 
of fluctuations, in accordance with a cosmological model on 
a uniform grid which covers the whole computational box. 
The fluctuation field is perturbed in Fourier space and the 
Zel'dovich displacement is assigned to the particles. The grid 
defines the lowest level of resolution of the simulation. Subse- 
quently, simulations are run with a fixed number of particles 
and halos arc identified in the sirrmlatioii. Once marked, the 
selected halos are re-simulated with higher mass and force 
resolution. 

Here, we propose simulations with a different approach. 
Assuming a spherical symmetry, the dark matter initial dis- 
tribution at z=100 of a single halo is very close to a uniform 
density distribution in expansion. A small negative radial 
density gradient causes each shell to reach the maximum 
expansion at the dynamical time corresponding to its inner 
average density. Because of the negative density gradient, 
the dynamical time is higher in the most external shells. 
This fact illustrates why the initial negative density gradi- 
ent is equivalent to a mass aggregation history (MAH) for 
the halo. In our case, for a given halo present mass, a MAH is 
obtained from the EPSF developed by Avila-Reese, Firmani 
and Hernandez (1998; hereafter AFH). Using the dynami- 
cal time of each shell implicit in the MAH we recover the 
initial density gradient, i.e. the initial condition at z=100. 
The stochastic nature of the density fluctuations induces a 
stochastic structure on the MAH. We have limited our anal- 
ysis on the average MAH from which a given present halo 
mass has grown. This initial condition for the description 
of the cosmological evolution of growing halos offers clear 
advantages with respect to the simulations of isolated halos 
starting with a NFW or a Hearnquist density profile. How- 
ever, because the mass accretion is described in a spherical 
symmetric frame, in this scheme the tidal destructive action 
of strong mergers on a central dense halo core is not taken 
into account. This fact establishes a limit on our approach 
that works in the sense to make more difficult the forma- 
tion of inner shallow density profiles. Because of our goal is 
to explore the capability of self-interaction to create shal- 
low cores, our results should be strengthened if a detailed 
cosmological simulation is made. Our result will be partic- 
ularly correct in the case of dwarf and LSB galaxies which, 
being fragile disks and field objects, their MAH could not 
be affected by major mergers. In the case of galaxy clusters 
our result is purely indicative and obeys to a completeness 
criterion. In fact in this case major mergers may amplify 
the formation of shallow cores as well as they may introduce 
significant deviations from the spherical symmetry. Taking 
into account the previous discussion and considering that 
our goal is to analyze the ability of self-interaction to create 
shallow cores, we concentrate our attention on the average 
MAH. The lack of any stochastical information from ob- 
servations makes unnecessary any analysis which takes into 
account the stochastical nature of MAHs. 

The following steps are taken: 

• We Eissign to the halo an initial density profile described 
by a collection of concentric shells. For each shell, the MAH 
of the halo gives: the time of maximum expansion tmax, the 



radius of the maximum expansion, the initial radius at a 
given rcdshift, and the cumulative mass within the shell. 

• The dynamical evolution of the halo within a non- linear 
regime is followed by the N-body code. 

At the beginning of our N-body simulations the radius of 

each particle is interpolated to the initial radius, determined 
by the halo MAH at a cumulative value of the mass. We have 
used different techniques to distribute the particles (the halo 
mass) within the shells at the beginning of the simulation, 
with a grid or with isotropic random direction. In the first 
case, particles are assigned on a spherical mesh points and 
are originally equally spaced on the mesh, then the distance 
of each particle from the centre is computed according to 
the halo MAH. In the second case, for each particle the rar 
dius is isotropically assigned using a Monte Carlo routine, 
while the magnitude of the radius is fixed by the cosmo- 
logical MAH. Both methods in assigning the particles at 
the initial conditions have limits and advantages. The first 
method uses a grid to assign particles and thus it is closer to 
the standard assignment used in the original version of the 
code. We have used it for the cosmological collisionless and 
coUisional runs. However, in some cases, a weak filament is 
seen in the virialized halo. This elongated structure is an 
artificial numerical effect, and some tests have shown that 
it is present already in the original version of the code. On 
the other hand, the second technique prevents the formation 
of elongated structures but since the particles are assigned 
with random isotropic directions, they rapidly clump, due to 
the gravity, to form substructures at very early times that 
survive during the virialization of the halo. These substruc- 
tures are spurious. Thus, for the cosmological simulations 
we prefer to assign particles to a grid. The Hubble fiow at 
z = 100 was assigned to the particles, as initial velocity. The 
initial value of z has been chosen with the criterion to avoid 
any spurious formation of soft cores. 



3 MODELLING COLLISIONLESS CDM HALOS 

The focus of the approach proposed here is on the explo- 
ration of density profiles of halos of any size and mass. We 
show now that by making use of the N-body code with 
the initial conditions assigned with our new numerical tech- 
nique, we can reproduce the NFW density profile for a single 
collisionless dark matter cosmological halo of any mass, af- 
ter the code is run for a few hours (CPU time) on a standard 
workstation. 

We carried out simulations of galactic halos as well 

as galaxy cluster halos within a fiat ACDM universe with 
matter density fim = 0.3, cosmological constant JIa = 0.7, 
and rms linear fluctuation amplitude in Sh~^ Mpc spheres 
of f78 = 1. Each simulation followed the trajectories of 
N = 80, 000 dark matter particles within a physical cube 
of side 0.5 Mpc from z=100 to the present. The Hubble 
constant is assumed to be of 65 km s^^ Mpc^^ (fe=0.65). 
Tests with more particles did not change the results. One re- 
finement level was introduced in the region of highest mass 
density. 

On galactic scales, for a halo mass of IQ^^h^^ Mq the 
mass per particle used is = M/80, 000 = 1.25 • l(fh~'^ 
Mq and the softening length is fixed in order to obtain a 
spatial effective resolution of « 1 kpc. On scales of 
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galaxy clusters we have carried out simulations of a sin- 
gle halo with M = lO'^^ft"^ Mq with a mass resolution of 
rup = 1.25 • 10^°/i~^ Mq and an effective spatial resolution 
of « 20 fe-i kpc. 

Wc stress that the cluster simulation was obtained in 
the same cube of side 0.5 Mpc. In fact, as gravitation 
is a scale invariant physical process, a halo of amy mass and 
size may be simulated in the same fixed volume. Of course, 
because of the hierarchical process of mass accretion, the 
mass aggregation history of a galactic halo is different from 
the galax;y cluster halo merging history. Thus, once the MAH 
is generated for a halo of galaxy cluster size, we can evolve 
the dark halo by N-body simulations by rescaling the mass 
M and the radius R to any mass Mq and initial radius Ro 
with the following relations: 



M = mMo, 

R = rRo, 

t = ^o(rVm)°■^ 



v'^ = vo {m/r) , 



(1) 
(2) 
(3) 
(4) 



where m, r are scaling parameters. This work simulates 

galaxy cluster halos by rescaling the cluster mass to mass 
Mo of a galactic halo with lO" /i^^M©. 

Simulated density profiles of coUisionless dark matter 
are shown in the top panel of Figure 1 for two halos with 
the following masses: M = 10" H-^Mq (filled circles) and 
M = 10^^ h'^Mpj (open circles). Solid lines are the NFW 
profiles for the same halos. The good agreement between 
the halo density profiles obtained with our method and the 
NFW model is encouraging. The velocity dispersion radial 
profiles corresponding to the same halos are drawn in the 
bottom panel of Figure 1 with filled symbols corresponding 
to the galactic halo of M = 10^^ h~^MQ and open symbols 
to galaxy cluster halo of M = 10^^ h~^MQ, as derived by 
our simulations. The overlapped solid lines are the radial 
velocity dispersion for the NFW model with the same mass. 

A more quantitative comparison to the NFW model for 
the two halos is obtained deriving values for the concen- 
tration parameter, defined as the ratio between the virial 
radius and the scale radius c = r^ir/rs (Navarro, Frenk 
& White 1997). For the halos of M = lO" H'^Mq and 
M = 10^^ Mq we estimate c=9.5 and 5, respectively 
for our realizations. The concentration values of our run are 
plotted in the Figure 2 as filled points and compared to the 
mass-concentration relation of the NFW model (dashed line) 
*. The halos are a little more concentrated than the aver- 
age concentration predicted by the NFW model. However, 
high-resolution N-body simulations have shown a spread 
around the average mass-concentration relation (Wechsler 
et al. 2002), thus the concctrations we find in our two real- 
izations are inside the spread. 



* The mass-concentration relation of the NFW model is derived 
using the program kindly made available to the community by 
J.Navarro. 
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Figure 1. In the top panel simulated density profiles of coUi- 
sionless dark matter are shown using non-standard N-body sim- 
ulations for halos with masses: M= 10^^h~^ Mq (filled circles) 

and M= lO^^h^^ AIq (open circles). A Hubble constant of 65 km 
s^-'^ Mpc^-*^ is adopted. Solid lines are the NFW profiles for the 
same halos. Corresponding velocity dispersion radial profiles are 
drawn in the bottom panel with filled symbols corresponding to 
the galactic halo and open symbols to the galaxy cluster halo. 
The solid lines are the radial velocity dispersion for the NFW 
model. 



4 MODELLING COLLISIONAL CDM HALOS 

The idea of a collisional dark matter suggested by Spergel 
& Steinhardt (2000), is capable of producing soft cores and, 
at the same time, preserves the success of CDM models in 
explaining the observed properties of the universe at large 
scales. In fact, in the hierarchical framework of structure 
formation, the halo continuously accretes mass by merging. 
However, collisions between particles cause heat transfer in- 
wards. Since its negative heat capacity, the core gains energy 
and expands, producing a soft core. However, the idea of 
weak self-interacting dark matter poses the following ques- 
tion. What is the final equilibrium configuration of the halo? 
The initial state is a halo described by a NFW model far 
from thermal equilibrium, and the final state is a halo with 
a thermalized core. 

We explore the dynamical properties of collisional dark 
matter studying isolated or cosmological halos of any mass 
using N-body techniques. We use a Monte Carlo technique 
to implement the self-interaction of the dark matter par- 
ticles in coUisionless N-body code. The scattering process 
between the dark particles is implemented with the follow- 
ing algorithm. Within a timestep At, the probability Pi of 
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Figure 2. Concentration parameters c versus the halo mass for 
the two halos of M= IQ^^H-^Mq and M= IO^^H-^Mq as ob- 
tained from our runs (filled points). The dashed line is the halo 
Mass-concentration relation of the NFW model. 



a particle with velocity Vi to interact with another particle 
is: 

Pi = Ti At, (5) 
with Ti the scattering rate for a particle with velocity Vi: 



ri = s 



AT / cr \mj [Vj — Vi\ 



V 



(6) 



where Tifnij/V is the estimate of the local density, with 
V = 47rr^/3. Although a constant cross section is used by 
several authors, here we assume a cross section inversely 
proportional to the particle velocity dispersion as found to 
be justified by the observations in the central densities of 
gala^xies and galaxy clusters (Firmani et al. 2001; Wyithe et 
al. 2001) 



■ V ~ const 



(7) 



Within a timostop the algorithm is as follows: 



i) For each particle the interaction probability Pi is com- 
puted. 

ii) For each scattering event the collision is with one of the 
nearest particles (N=16, but increasing N does not change 

the result). 

in) Once a particle is selected for the interaction, the scat- 
tering is isotropic and the diffusion is computed. 

Concerning the diffusion process, for particles with centre- 
of-mass velocity vh = {vl + '>^j)/2 and velocity difference 



Figure 3. In the top panel the time evolution of the density 
profile is shown for a halo of mass M= 10"/i-1Mq described 
by a singular isothermal sphere. This run uses coUisionless dark 
matter particles. After several timeseales, the density distribution 
is represented by the solid line, while the short dashed line is 
the profile at the beginning of the simulation (overlapped). In 
the bottom panel the corresponding radial velocity dispersion is 
shown. 



= {vi — Vj)/2, a random direction e is selected and new 
velocities are assigned to the colliding particles i and j 



Vi=VB + -r-flM 
e 



Vj =vb- tt ^0 . 



(8) 



(9) 



The energy and linear momentum are preserved. The nu- 
merical method is similar to the algorithm implemented by 
Burkert (2000) and Yoshida et al. (2000) with the differ- 
ence that the cross section is dependent upon the particle 
velocity. For the cross section we assume: 



/ o- \ .^^-24/ 100 km s-^ \ cm^ 
\mxJ \ V J 



GeV 



(10) 



In order to simulate a galaxy cluster halo with coUisional 
dark matter, the cross section value was rescaled as well as 
the radius, velocity, time, and mass according to eqs.(4) to 
(6) inclusive. 



5 ISOLATED HALOS AND CORE COLLAPSE 

For isolated halos the balance between continuous mass ag- 
gregation due to the halo merging history and the thermal- 
ization process induced by collisions fails. In this section we 
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investigate wliether tliis case can produce core collapse with 
a consequent gravothermal catastrophe in a short time. 

First, we test if the code is able to preserve the state 
of dynamical equilibrium in a model. We have simulated a 
halo with colUsionless dark particles and a mass profile of 
an isothermal sphere. This choice is related to the known 
properties of this model. In the top panel of Figure 3, the 
time evolution of the density profile is shown for several dy- 
namical timescales. The density profile preserves the shape 
for all times. After several timescales, the density distribu- 
tion is represented by the solid line, while the short dashed 
line is the profile at the beginning of the simulation. In the 
bottom panel the corresponding radial velocity dispersion is 
shown. 

Note that the isothermal sphere was settled by assigning 
to the particle velocities with amplitudes consistent with the 
Maxwellian velocity distribution function, f{v) of the model. 
In fact, in this case the initial profile is at equilibrium and 
does not show any adjustment during dynamical evolution. 
The result is remarkable: the dynamical equilibrium of the 
model is preserved for all times. The only modification to the 
profile is after many dynamical timescales in the outer parts 
due to a evaporation process of the particles. Furthermore, 
this model tests the softening radius fixed in our simulations. 
For a galactic halo of 10^^ h~^MQ the softening radius was 
fixed in order to obtain a spatial effective resolution of 1 
kpc. 

We have followed the evolution of dark density profiles for 
a galactic halo of 10^^ h~^MQ formed by collisional dark 
matter and characterized by the following initial mass dis- 
tribution: the Hernquist model and the King profile. 



5.1 The Hernquist profile 

An isolated galactic halo model described by an initial Hern- 
quist profile (Hernquist 1990) has attracted considerable at- 
tention in the study of self-interacting dark matter. 

For this model characterized by a central density 
cusp (1/r) similar to a NFW profile, Burkert (2000) and 
Kochanek & White (2000) axgue that soft cores are pro- 
duced in a short dynamical timescale induced by collisions. 
However the shallow cores so created are unstable and a sin- 
gular isothermal sphere is produced as a final halo density 
profile. 

We set up the initial conditions as a Hernquist pro- 
file with the halo mass of lO'^'^ h~^MQ and a characteristic 
scale radius a = kpc. Following the analytical pro- 

file, particle positions are distributed isotropically and the 
magnitude of the velocity is computed from the energy dis- 
tribution function, f{E), of this model, also assigned with 
an isotropic distribution. 

First, we run a model without collisions between parti- 
cles and we have shown that the shape of the profile is not 
modified during dynamical evolution. The dynamical equi- 
librium is preserved. The dynamical time is defined: 



tdy 



( — )■ 



(11) 
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thus, the dynamical time is = 6.7 Gyr (t^n is referred 
to the Hernquist model). In the top- left panel of Figure 4, 
the time evolution for the Hernquist profile is shown for self- 



Figure 4. Left panels. The time evolution for the Hernquist pro- 
file is shown for self-interacting cross section value obtained as- 
suming a cross section inversely proportional to the halo disper- 
sion velocity, {(j/mx) ■ vioo = 10"^'' cm'^ GcV"-', where dioo 
is the halo dispersion velocity in units of 100 km s^^ . The ini- 
tial Hernquist profile is shown as a short-dashed line. After 2.25 
dynamical timescales a minimum central density is reached (long- 
dashed lines) and the core collapse after 5 dynamical timescales 
(solid lines). Corresponding velocity dispersion radial profiles are 
shown in the bottom. Right panels. The time evolution of the 
King profile is plotted. Dashed lines are the initial King profile 
(t=0). After a while, the initial flat core reduces its size towards 
core collapse. In the plot the time evolution halo density profile 
is shown, from the lower solid line to the higher solid line rep- 
resenting the profile after 0.8 t^.^- In the bottom panel of the 
same figure the radial velocity dispersion profile is shown for the 
model. 



interacting cross section value obtained assuming a cross 
section inversely proportional to the halo dispersion veloc- 
ity, {a/rrix) ■ I'loo ~ 10"^'' cm^ GeV""*^, where «ioo is the 
halo dispersion velocity in units of 100 km s^^. The ini- 
tial Hernquist profile is shown as a short-dashed line. After 
2.25 dynamical timescales a minimum central density value 
is reached (long-dashed line), increasing after 4 dynamical 
timescales and reaching a lack of the core after more than 
5 dynamical timescales (solid line) , due to the collisions be- 
tween particles. We also run simulations with very high cross 
section values, finding that once the soft core is formed, it 
disappears in a short time, in agreement with the work of 
Kochancck & White (2000) and Burkert (2000), when the 
cross section was adopted to be independent of the parti- 
cle collision velocity. Thus, high cross section values seem to 
work on a catalyst of the core collapse. 
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5.2 The King profile 

We have simulated an isolated halo of 10^^ Mq with a 
King model mass distribution. In fact, the dark density pro- 
file of LSB and dwarf galaxies is well matched by a King 
profile with form parameter P = 8 (Firmani et al. 2001). 
Wo test whether starting with a isothermal configuration 
and a shallow core the isolated halo undergoes core collapse 
due to the presence of collisions between the dark particles. 
We want to follow the halo core evolution through isother- 
mal equilibrium states. After assigning isotropic positions 
to the particles, wc pay attention that at the beginning of 
the simulation the profile is at equilibrium, and we calculate 
the velocity magnitude from the King distribution function 
f{v). Again, particle velocities are isotropically distributed 
with a Monte Carlo routine. 

First, we have run the model with a/rux = 0. The shape 
of the King profile is preserved at all times. We now run the 
model in a strong self- interacting regime: {a/rrix) ■ uioo — 
5- 10^'^^ cva? GeV^^. In fact, for weak cross sections, we have 
verified that this model evolves slowly towards core collapse. 
In this large cross section value is assumed in order 

to hasten core collapse. The dynamical time depends upon 
the adopted concentration for the King model. Defining the 
dynamical time as above, we find that t^^^ ~ 4 t^y^ (where 

is referred to the King model and t^y^ to the Hernquist 
model, respectively). In the top-right panel of Figure 4, the 
time evolution of the King profile is plotted. Dashed lines are 
the initial King profile (t=0). After a while, the initial flat 
core reduces its size towards core collapse. In the plot the 
time evolution halo density profile is shown, from the lower 
solid line to the higher solid line representing the profile after 
0.8 td^n- In the bottom panel of the same figure the radial 
velocity dispersion profile is shown for the model. 

It is interesting to note that any reasonable value for 
the cross section (weak or strong) produces core collapse. 
Weaker coUisional regimes will produce core collapse in a 
longer time. Of interest is the evidence that a dark halo 
described by a King profile evolves via equilibrium states 
towards the gravothermal catastrophe under the effect of 
collisions. 



6 COLLISIONAL HALOS IN A 

HIERARCHICAL COSMOLOGICAL 
FRAMEWORK 

The question we investigate in this section is: assuming a 
modest self-interacting cross section with {a jm^') ■ wioo = 
10"^"' cm^ GeV^^ are halo soft cores produced by coUisional 
dark matter in a cosmological framework stable or can they 
collapse? In other words, can dark halo density profiles be 
soft at the centre and changing slowly their slope reaching 
the NFW shape or steeper, due to the collision, in a Hubble 
time? 

Using the N-body technique described above we run 
simulations for a galactic halo and a cluster halo. In Fig- 
ure 5 we show the dark matter density profiles of a halo of 
M= 10" H'^Mq as obtained in a coUisional CDM model. 
The halo have been modelled assuming a energy dependent 
cross section with the same value as above. In the top panel 




Q.5 1 
LogR [kpc h~ 



Figure 5. In the tup panel the dark matter density profiles of a 
halo of M= lO^'^ h~^MQ are shown as obtained in a coUisional 
CDM scenario. The halo have been modelled assuming a energy 

dependent cross: (cr /mx)-vioQ = 10"^'* cm^ GcV^"'". Short dashed 
lines are the halo density profile after 7.5 Gyr, wfiilc long-dashed 
lines are the lialo mass density distribution after a Hubble time. 
In tiie same plot the solid line is the NFW profile. In the bottom 
panel the corresponding halo radial dispersion velocity profiles 
are shown. 



the short dashed line shows the halo density profile after 
7.5 Gyr, while the long-dashed line is the halo mass density 
distribution after a Hubble time. In the same plot the solid 
line is the NFW profile. 

Note that within a Hubble time a modest self-interaction 
cross section value is able to create soft cores. No core col- 
lapse occurs within a Hubble time. 

The bottom panel of Figure 4 shows the corresponding 
halo radial dispersion velocity profiles. The velocity distri- 
bution of the particles reaches roughly a constant velocity 
dispersion under the effects of collisions resulting in a central 
non-singular isothermal density profile. 

In Figure 6 the dark density profile is shown for a halo 
of 10^'' /i^^Mq. The long-dashed line is the halo density 
profile after one Hubble time and the solid line is the NFW 
profile. Even in this case the halo does not undergo the core 
catastrophe within a Hubble time. We interpret this result 
assuming that this is the case in which any trend towards the 
gravothermal catastrophe is avoided by the competition be- 
tween a mass aggregation determined by a continuous halo 
merging history and a thermalization process by collisions. 
In fact wc assign to the particles cosmological initial condi- 
tions assuming that the halo merging history is a continuous 
process. 

In Figure 7, the central densities and the core radii 
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Figure 6. The dark density profile is shown for a halo of 
10^^ h~^MQ, using a energy independent cross section: {c/mx) ■ 
vioo ~ 10~^* cm^ GeV~^. The long-dashed line is the halo den- 
sity profile after one Hubble time and the solid line is the NFW 
profile. 



predicted by the models are compared with those inferred 
from HI rotation curves of dwarf and LSB galaxies (diagonal 
crosses), Ha rotation curves from Swatcrs ct al. (2000), tak- 
ing into account the halo adiabatic contraction (open circles) 
and Ha rotation curves of LSB and dwarf galaxies (filled 
circles) by Marchesini et al. (2002). The filled circles on the 
right of the panels are two galaxy clusters with evidence 
of soft cores: CI 0024+1654 and Abell 1795. The dashed 
lines are the predictions of the model for a cross section 
(fx/mx) • wioo = lO"^"* cm^ GeV-^ 

A feature of self-interacting dark matter, with any cross 
section is to produce halos that arc more spherical than the 
CDM halos, due to the isotropic nature of collisions. Thus, 
collisions tend to isotropize the velocity ellipsoid. Conse- 
quently the shapes of self-interacting dark matter halos and 
CDM halos should be different and provide an interesting 
observational constraint to discriminate the two scenarios. 
Miralda-Escude (2002), providing lensing data for the galaxy 
cluster MS2137-23, shown that at a radius Ji « 70 kpc, 
anisotropies start in the velocity ellipsoid, constraining the 
cross section between dark particles to be less than 0.02 
cm^g~^ in order to create anisotropies in agreement with 
the observations. 

Dave et al. (2001) studied the the axisymmetry and 
the flattening for galactic halos using N-body simulations 
with constant cross section between dark particles. After the 
completition of this work, a study of N-body simulations of 
self-interacting dark halos appeared by Colin et al. 2002, in 



Figure 7. The halo central densities and the core radii predicted 
by the models are compared with those inferred from HI rotation 
curves of dwarf and LSB galaxies (diagonal crosses) , Ha rotation 

curves from Swatcrs ct al. (2000), taking into account the halo 
adiabatic contraction (open circles) and Ha rotation curves of 
LSB and dwarf galaxies (filled circles) by Marchesini et al. (2002). 
The filled circles on the right of the panels arc two galaxy clusters 
with evidence of soft cores: CI 0024+16.5-1 and Abcll 1795. Dashed 
lines are the halo central densities and core radii as a function of 
the maximum circular velocity predicted by self-intera<;ting dark 
matter with {ajmx) ■ I'loo — 10~^* cm^ GeV"'^. 

which the authors analyse the ellipticities in the halos, even 
using a cross section inversely proportional to the halo dis- 
persion velocity. Due to the limitations of our method, cos- 
mological triaxiality and a comparison to the lensing data 
cannot be developed within our approach. Hovowcr the ef- 
fects of the self-interaction on an elongated structure artifi- 
cially introduced in the early evolution of the halo have been 
explored. We define the axisymmetry q and the flattening s, 
following the same prescription used in Dave et al. (2001). 
A inertia tensor is defined: 




with the sum is over all particles with coordinates 
(a;i, a:2, 0:3) and distance r from the halo centre and s and q 
the axis ratio: 



with Mxx^ MyyjM^z are the eingenvalues of the inertia ten- 
sor M. 



© 2002 RAS, MNRAS 000, 1-10 



The Halo Density Profiles with Non-StandardN-hody simulations 9 



1.2 



"T — \ — I — I — I — r 



M=10'5 




J \ I I I I \ I L 



-2 -1.5 -1 -0.5 
Log r/r , 



Figure 8. Axisymmetry q and flattening s for a cluster lialo 
of M=10^^ h~^MQ versus r/ruir- Solid lines represent s (lower 
line) and q (higher line) runs with collisionless dark matter halos. 

Dashed lines are simulations with coUisional dark matter with 
cross section {a /rax) ■ «^ioo — 10^"^'* cm-^ GeV~^. 



Hernquist model and the King profile. Surprisingly, the viri- 
alized halo density profiles for all models undergo core col- 
lapse induced by collisional dark matter. Once the core size 
is reduced, core collapse is a very rapid process. The shape 
of the cross section (constant or dependent upon the particle 
collision velocity) does not affect core collapse. Our N-body 
simulations performed with a cross section dependent on the 
particle velocity dispersion produce similar results to the N- 
body simulations by Burkert (2000) and Kochanek & White 
(2000) obtained with a constant cross section. 

• In the hierarchical framework, the halo core catastro- 
phe is avoided from the balance between mass aggregation 
process induced by gravitation and thcrmalization process 
induced by collisions. If the balance fails, collisions between 
particles induce core collapse in a short time. Large cross 
section values cause even faster core collapse. 

• If the evidence of soft cores is accepted for CI 
0024-1-1654 and Abell 1795, the observed central densities of 
both clusters and dwarf galaxies are reproduced assuming 
a cross section value of: {a /nix) ■ fioo ~ 10~^* cm^ GeV""*^. 
The shape and the value of the cross section suggested by 
our work are in the range of values in which the core col- 
lapse and the evaporation problem are avoided (Gnedin & 
Ostrikcr 2001; Hcrmawi & Ostriker 2002). 

Observations of the mass distribution at the centre of galaxy 
clusters are crucial for the debate regarding the nature of 
the dark matter. If future observations cannot confirm the 
evidence for soft cores on galaxy cluster scales, then the soft 
core question remains confined to dwarf and LSB galaxies. In 
this case different mechanisms could play a role in producing 
soft cores in dark matter dominated galaxies. 



Since the observational constraints on anisotropies ex- 
istence come from galaxy cluster we feel more interesting to 
analyse the galaxy cluster scale halo with M=10^^ /i~^M© 
and to assume the same cross section capable to repro- 
duce the soft core in agreement with the observations of 
C10024-hl654, {o/mx) ■ wioo ^ 10"^^ cm^ GeV~\ Figure 
8 shows how the galaxy cluster halo is more spherical if 
a self-interacting scenario is working. For a cluster with 
Tvir = 1.5 Mpc any early triaxiality desappears inside a re- 
gion of roughly 120 kpc (for h=0.65). This is not surprising 
because particle collisions induced by self-interacting dark 
matter tend to recover the spherical symmetry. However in 
realistic cases major mergers may be a continuous source for 
anisotropies, even within a few hundred kpc from the centre. 



7 DISCUSSION AND CONCLUSION 

We have explored the dark halo density profiles in a colli- 
sional dark matter scenario studying isolated or cosmologi- 
cal halos of any size and mass. This study presents two key 
points: we use a modified N-body code where the cosmo- 
logical initial conditions are assigned with a new approach; 
and ii) wc use a self-interacting cross section inversely pro- 
portional to the particle collision velocity. 

Our results may be summarized as follows: 

• Concerning the dynamical study of isolated halos we 
have followed the evolution of different initial profiles: the 



ACKNOWLEDGMENTS 

We wish to thank the anonymous referee for his comments 
and suggestions that have improved this work. We are grate- 
ful to Hugh Couchman for having made available his adap- 
tive P^M-SPH code HYDRA. ED thanks Henry Lee for an 
early reading of the paper. 



REFERENCES 

Arabadjis, J.S., Bautz, M.W., Garmire, G.P. 2002, ApJ, 572, 66 
Avila-Reese, V., Firmani, 0., Hernandez, X. 1998, ApJ, 505, 37 
Burkert A., 1995, ApJ, 447, L25 
Burkert, A. 2000, ApJ, 534, L143 

Colin P., Avila-Reese V., Valenzuela O., Firmani C, 2002, ApJ, 

submitted, (astro- ph/0205322) 
Couchman, H.M.P., Thomas, P.A., Pearee, F.R. 1995, ApJ, 452, 

797 

Dave, R., Spergel, D.N., Steinhardt, P.J., Wandelt, B.D. 2001, 
ApJ, 547, 574 

David L.P., Nulsen P.E.J. , McNamara B.R., Forman W., Jones 
C., Ponman T., Robertson B., Wise M. 2001, ApJ, 557, 456 

de Blok, W.J.G., McGaugh, S.S., Rubin, V.C. 2001, AJ, 122, 2396 

Ettori, S., Fabian, A. 2002, Chemical Enrichment of Intraclus- 
ter and Intergalactic Medium, ASP Vol. 253, ed. by Fusco- 
Femiano, F. Matteucci, 2002., p.23 

Firmani, C, D'Oiigliia, E., Chincarini, G., Hernandez, X., Avila- 
Reese, V. 2001, MNRAS, 321, 713 

Flores, R.A., and Primack, J., 1994, ApJ, 427, LI 

Fukushige T., Makino J., 2001, ApJ, 557, 533 



© 2002 RAS, MNRAS 000, 1-10 



10 E. D'Onghia}''^, C. Firmani^'^, G. Chincarini^'^ 



Gnedin, O.Y., & Ostrikcr, J. P. 2001, ApJ, 561 
Hennawi, J.F., Ostriker, J. P. 2002, ApJ, 572, 41 
Hernquist, L. 1990, ApJ, 356, 359 

Hogan, C.J., & Dalcanton, J.J. 2000, Phys.Rev., D62, 063511 

Kochanek, C.S., White, M., 2000, ApJ, 543, 514 

Maxchesini, D., D'Onghia, E., Chincaxini G., Firmani, G., Gon- 

coni. P., Molinari, E., Zacchei, A. 2002, ApJ, 572, 66 
Miralda-Escude, J. 2002, ApJ, 564, 1019 
Moore B., 1994, Nature, 370, 629 

Moore B., Quinn T., Govcrnato F., Stadcl J., Lake G., 1999, 

MNRAS, 310, 1147 
Navarro J., Frenk G.S., White S.D.M., 1997, ApJ, 490, 493 
Ostriker, J.P. 2000, Phys.Rev.Lett., 84, 23 
Salucci, P., Burkert, A. 2000, ApJ, 537, L9 
Sellwood, J.A. 2000, ApJ, 540, LI 

Spergel, D.N., Steinhardt, P.J., 2000, Phys.Rew., 84, 3760 
Swaters, R.A., Madore, B.F., & Trewella, M. 2000, ApJ, 531, L107 
Tyson, J.A., Kochanski, G.P., dell' Antonio, LP., 1998, ApJ, 498, 

L107 

van den Bosch E.G. & Swaters R.A., 2001, MNRAS, 325, 1017 
Wechsler R.H., Bullock J.S., Primack J.L., Kravtsov A.V., Dekel 

A., 2002, ApJ, 568, 52 
White S.D.M. and Rees M. 1978, MNRAS, 183, 341 
Wyithe, J.S.B., Turner, E.L., Spergel, D.N. 2001, ApJ, 555, 504 
Yoshida, N., Springel, V., White, S.D.M., Tormen, G. 2000, ApJ, 

544, L87 

This paper has been produced using the Royal Astronomical 
Society/Blackwell Science style file. 



© 2002 RAS, MNRAS 000, 1-10 



